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Abstract. Large-scale controlled evacuations require emergency services to se¬ 
lect evacuation routes, decide departure times, and mobilize resources to issue 
orders, all under strict time constraints. Existing algorithms almost always al¬ 
low for preemptive evacuation schedules, which are less desirable in practice. 
This paper proposes, for the first time, a constraint-based scheduling model that 
optimizes the evacuation flow rate (number of vehicles sent at regular time in¬ 
tervals) and evacuation phasing of widely populated areas, while ensuring a non- 
preemptive evacuation for each residential zone. Two optimization objectives are 
considered: (1) to maximize the number of evacuees reaching safety and (2) to 
minimize the overall duration of the evacuation. Preliminary results on a set of 
real-world instances show that the approach can produce, within a few seconds, a 
non-preemptive evacuation schedule which is either optimal or at most 6% away 
of the optimal preemptive solution. 

Keywords: constraint-based evacuation scheduling - non-preemptive schedul¬ 
ing - phased evacuation - simultaneous evacuation - actionable plan - real-world 
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1 Introduction 

Evacuation planning is a critical part of the preparation and response to natural and 
man-made disasters. Evacuation planning assists evacuation agencies in mitigating the 
negative effects of a disaster, such as loss or harm to life, by providing them guidelines 
and operational evacuation procedures so that they can make informed decisions about 
whether, how and when to evacuate residents. In the case of controlled evacuations, 
evacuation agencies instruct each endangered resident to follow a specific evacuation 
route at a given departure time. To communicate this information in a timely fashion, 
evacuation planners must design plans which take into account operational constraints 
arising in actual evacuations. In particular, two critical challenges are the deployment of 
enough resources to give precise and timely evacuation instructions to the endangered 
population and the compliance of the endangered population to the evacuation orders. In 
practice, the control of an evacuation is achieved through a mobilization process, during 
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which mobilized resources are sent to each residential area in order to give instructions 
to endangered people. The number of mobilized resources determines the overall rate at 
which evacuees leave. Finally, to maximize the chances of compliance and success of 
a controlled evacuation, the evacuation and mobilization plans must be easy to deploy 
for evacuation agencies and should not leave, to the evacuees, uncontrolled alternative 
routes that would affect the evacuation negatively. 

Surprisingly, local authorities still primarly rely on expert knowledge and simple 
heuristics to design and execute evacuation plans, and rarely integrate human behavioral 
models in the process. This is partly explained by the limited availability of approaches 
producing evacuation plans that follow the current practice. Apart from a few exceptions 
II2I5I6I9I111131 existing evacuation approaches rely on free-flow models which assume 
that evacuees can be dynamically routed in the transportation network 03I1OI14I . These 
free-flow models however violate a desirable operational constraint in actual evacuation 
plans, i.e., the fact that all evacuees in a given residential zone should preferably follow 
the same evacuation route. 

Recently, a handful of studies considered evacuation plans where each residential 
area is assigned a single evacuation path. These studies define both a set of evacuation 
routes and a departure schedule. Huibregtse et al. 0 propose a two-stage algorithm that 
first generates a set of evacuation routes and feasible departure times, and then assigns 
a route and time to each evacuated area using an ant colony optimization algorithm. In 
subsequent work, the authors studied the robustness of the produced solution ||7], and 
strategies to improve the compliance of evacuees ISl. Pillac et al. ifTJI first introduced 
the Conflict-based Path Generation (CPG) approach which was extended to contraflows 
by Even et al. 0. CPG features a master problem which uses paths for each residential 
node to schedule the evacuation and a pricing problem which heuristically generates 
new paths addressing the conflicts in the evacuation schedule. 

These evacuation algorithms however do not guarantee that evacuees will follow 
instructions. If the evacuation plan contains forks in the road, evacuees may decide to 
change their evacuation routes as the evacuation progresses. This issue is addressed in 
ifTTel which propose evacuation plans without forks. The resulting evacuation plan can 
be thought of as a forest of evacuation trees where each tree is rooted at a safe node (rep¬ 
resenting, say, an evacuation center) and with residential areas at the leaves. By closing 
roads or controlling intersections, these evacuation trees ensure the compliance of the 
evacuees and avoid congestions induced by drivers slowing down at a fork. Even et al. 
0 produce such convergent evacuation plans by decomposing the evacuation problem 
in a tree-design problem and an evacuation scheduling problem. Andreas and Smith m 
developed a Benders decomposition algorithm that selects convergent evacuation routes 
that are robust to a set of disaster scenarios. 

All of the approaches reviewed above allow preemption; The evacuation of a resi¬ 
dential area can be interrupted and restarted arbitrarily. This is not desirable in practice, 
since such schedules will confuse both evacuees and emergency services, and will be 
hard to enforce. Non-preemptive schedules have been considered in 141121 in rather 
different ways. In ||4l, a phased evacuation plan evacuates each area separately, guar¬ 
anteeing that no vehicles from different areas travel on a same path at the same time. 
By definition, phased evacuation does not merge evacuation flows, which is motivated 








by empirical evidence that such merging can reduce the road network capacities. The 
algorithm in ina is a column-generation approach for simultaneous evacuation, i.e., 
evacuations where multiple paths can share the same road segment at the same time. 
Each column represents the combination of a path, a departure time, and a response 
curve capturing the behavioral response of evacuees for each evacuation area. Here 
the flow rate of each evacuation area is restricted to pre-existing response curves, and 
columns are generated individually for each evacuation area. This column-generation 
approach requires a discretization of the evacuation horizon. 

This paper proposes, for the first time, a constraint programming approach to gen¬ 
erate non-preemptive evacuation schedules. It takes as input a set of evacuation routes, 
which are either chosen by emergency services or computed by an external algorithm. 
The constraint-based scheduling model associates a task with each residential area, uses 
decision variables for modeling the number of evacuees (i.e., the^ow), the number of 
vehicles to be evacuated per time unit (i.e., the ^ow rate), and the starting time of the 
area evacuation; It also uses cumulative constraints to model the road capacities. In 
addition, the paper presents a decomposition scheme and dominance relationships that 
decrease the computational complexity by exploiting the problem structure. Contrary 
to ifT^ . the constraint-programming model uses a decision variable for the flow rate of 
each evacuation area (instead of a fixed set of values) and avoids discretizing time. In 
contrast to the constraint-programming model allows for simultaneous evacuation 
while satisfying practice-related constraints. 

The constraint-programming model was applied on a real-life evacuation case study 
for the Hawkesbury-Nepean region in New South Wales, Australia. This region is a 
massive flood plain protected from a catchment area (the blue mountains) by the Warra- 
gamba dam. A major spill would create damages that may reach billions of dollars and 
require the evacuation of about 80,000 people. Preliminary experimental results indi¬ 
cate that the constraint-programming model can be used to generate non-preemptive 
schedules that are almost always within 5% of the optimal preemptive schedules gener¬ 
ated in prior work. These results hold both for maximizing the number of evacuees for 
a given time horizon and for minimizing the clearance time (i.e., the earliest time when 
everyone is evacuated). These results are particularly interesting, given that the optimal 
preemptive solutions produce evacuation plans which are far from practical. Indeed, 
Fig.[2shows the repartition of departure times for seven residential areas in the original 
HN80 instance using a preemptive schedule produced by the algorithm in ||6l. Observe 
how departure times are widely distributed within the scheduling horizon, indicating 
that the plan makes heavy use of preemption and is virtually impossible to implement 
in practice. Finally, experimental results on the phased version of the algorithm indicate 
that phased evacuations are much less time effective, and should not be the preferred 
method for short-notice or no-notice evacuation. 

The rest of this paper is organized as follows. Section [^defines the problem. Sec¬ 
tion [^presents the constraint-programming model, including the core model, the de¬ 
composition scheme, the dominance relationships, and the search procedure. Sectionj^ 
presents the experimental results. Section|^concludes the paper. 



Fig. 1; Departure times and flows of 7 residential areas of the HN80 instance with the 
preemptive algorithm FSP from Q. 

2 Problem Description 

The Evacuation Planning Problem (EPP) was introduced by the authors in ifTTll . It is 
defined on a directed graph Q = {Af = £ UT U S, A), where £, T, and S are the set 
of evacuation, transit, and safe nodes respectively, and A is the set of edges. The EPP is 
designed to respond to a disaster scenario, such as a flood, which may determine a time 
at which some edges become unavailable. Each evacuation node fc € £ is characterized 
by a number of evacuees dk, while each arc e is associated with a triple (te, Ug, be), 
where te is the travel time, rtg is the capacity, and he is the time at which the arc be¬ 
comes unavailable. We denote by e.tail (resp. e.head) the tail (resp. head) of an edge 
e. The problem is defined over a scheduling horizon H, which depends on the disaster 
forecast and the time to mobilize resources. The objective is either (1) to maximize the 
total number of evacuees reaching a safe node (for a fixed horizon) or (2) to minimize 
the time at which the last evacuee reaches a safe zone (for a variable horizon). In the 
following, we assume that the evacuation is carried out using private vehicles, but the 
proposed approach could be adapted to other contexts, such as building evacuation. 

This paper extends the EPP to the non-preemptive (simultaneous) evacuation plan¬ 
ning problem (NEPP) and the non-preemptive phased evacuation planning problem 
(NPEPP). Both are designed to assist evacuation planners with the scheduling of fully 
controlled evacuations. Given a set of evacuation paths, the NEPP decides the start 
time, flow, and flow rate at which vehicles are evacuating each individual evacuation 
node, ensuring that the evacuation operates without interruption. The NEPP allows sev¬ 
eral evacuation nodes to use the same road segments at the same time. In contrast, 
the NPEPP guarantees that no two evacuation nodes use the same road segment at the 
same time. The practical interest of the NPEPP is to evacuate designated priority areas 
quickly and efficiently, by limiting the risk of any delay caused by slowdown or traffic 
accidents which may result from merging traffic. 

Eormally, an evacuation plan associates with each evacuation area k £ £ exactly 
one evacuation path pk which is used to route all residents in fc to a same safe node. 
Let f2p^ = UfegfPfe *^he set of evacuation paths for all evacuations nodes in £. The 
characteristics of a path pk are as follows. We denote by Ap^ (resp. AfpA the set of 
edges (resp. nodes) of pk and by £{e) the set of evacuation areas whose path contains 
























edge e, i.e., e G Ap^. The travel time tk,n between the evacuation area k and a node 
n € Ku is equal to the sum of the path edges travel times separating k from n and 
tp^ is the total travel time between the start and end of p}~. The path capacity Up^ is the 
minimum edge capacity of pk- The last possible departure time LASTDEP(pfc) along 
path pk, i.e., the latest time at which a vehicle can depart on pk without being blocked, 
is easily derived from all tk,n{n G J^pk) the time he at which each path edge 
e € Ap^. becomes unavailable. If none of the path edges e G Apf. are cut by the disaster 
then LASTDEP(pfc) = oo; otherwise LASTDEP(pfe) = mineg^^^ {be — tk^e.head)- Note 
that the latest path departure time only depends on the path characteristics and not on 
H. 

3 The Constraint-Programming Model 

The NEPP and NPEPP are two optimization problems whose objective is either to max¬ 
imize the number of evacuees or to minimize the overall evacuation clearance time. The 
key contribution of this paper is to model them as constraint-based scheduling problems 
and to use CP to solve them. This modeling avoids time dicretization and makes it pos¬ 
sible to design non-preemptive plans with variable flow rates. This section presents the 
constraint-based scheduling models, including their decision variables, their domains, 
and the constraints common to both problems. This section then presents the constraint- 
based scheduling models for the NEPP in Sect. ?? and for the NPEPP in Sect. ??. 


3.1 Decision Variables 

The models associate with each evacuation area k G £ the following decision variables: 
the total flow of vehicles evacuated FLOWfc (i.e., the number of vehicles evacuated from 
area k), the flow rate ApLowk representing the number of vehicles departing per unit of 
time, the evacuation start time START^ (i.e., the time at which the first vehicle is evacu¬ 
ated from area k), the evacuation end time END^., and the total evacuation duration time 
DURfc. The last three decision variables are encapsulated into a task variable TASK^ 
which links the evacuation start time, the evacuation end time, and the evacuation dura¬ 
tion and ensures that STARTfc 4- DUR^ = END^. 

The decision variables range over natural numbers. The flow and flow rates can only 
be non-negative and integral since a number of vehicles is a whole entity. The models 
use a time step of one minute for flow rates and task variables which, from an opera¬ 
tional standpoint, is a very fine level of granularity: Any time step of finer granularity 
would only be too complex to handle in practice. The domains of the decision variables 
are defined as follows. The flow variable is at most equal to the evacuation demand: 
FLOWfc G [0, dfc] where [a, 6] = {u G N | a < v < b}. The flow-rate variable has 
an upper bound which is the minimum of the evacuation demand and the path capac¬ 
ity rounded down to the nearest integer, i.e., Aflow;, G [1j iiiin((ife, [up^J)]. The upper 
bounds for the evacuation start time and evacuation end time are the smallest of the 
scheduling horizon minus the path travel time, which is rounded up to the nearest inte¬ 
ger, and the latest path departure time, i.e., START^ G [0, min('H—[fpfc] , [LASTDEP(pfc)J )]. 
The evacuation of an area k can last at most dk minutes assuming the flow rate is set to 



one vehicle per minute; DUR^ G [0, dfc]. Note that the lower bound for duration is zero 
in order to capture the possibility of not evacuating the area. 

3.2 Constraints 

The NEPP requires to schedule the flow of evacuees coming from each evacuation 
area k on their respective path pk such that, at any instant t, the flow sent on all paths 
through the network does not exceed the network edges capacities. These flow con¬ 
straints can be expressed in terms of cumulative constraints. Consider an edge e and 
the set £ (e) of evacuation areas whose evacuation paths use e. For each evacuation area 
k G £(e), the model introduces a new task TASK| which is a view over task TASK^ 
satisfying: 

START^ = STARTfc -f ffc,e.ta*Z , DUR^ = DURfc , END^ = ENDfe -f ■ 

This new task variable accounts for the number of vehicles from evacuation area k trav¬ 
eling on edge e at any time during the scheduling horizon. Note that tk^e.taii is com¬ 
puted as the sum of the travel times on each edge, each rounded up to the next integer 
for consistency with the domains of the decision variables. While this approximation 
may slightly overestimates travel times, it also counterbalances possible slowdowns in 
real-life traffic, which are not taken into account in this model. 

The constraint-based scheduling model for the NEPP introduces the following cu¬ 
mulative constraint for edge e: 

cumulative({(TASK^, ApLowfc) | k G £(e)},Ue). 

The constraint-based scheduling model for the NPEPP introduces a disjunctive con¬ 
straint for edge e instead; 

dis junctive({TASK^ I fc € £’(e)}). (1) 

3.3 The Constraint-Based Scheduling Models 

We are now in a position to present a constraint-based scheduling model for NEPP-MF: 


max OBJECTIVE = ^ FLOWfc 


(2) 

k^S 



S.t. FLOW^® = DURfc X ApLoWfc 

VkG£ 

(3) 

FLOWfc = min(FLOWfc®, dfc) 

VkG£ 

(4) 

cumulative({(TASK^., ApLowfc) I k G £(e)},Ue) 

VeGA 

(5) 


The objective Q maximizes the number of evacuated vehicles. Constraints Q and Q 
link the flow, flow rate, and evacuation duration together, by ensuring that the total flow 
for each area k is the minimum of the evacuation demand and the flow rate multiplied 
by the evacuation duration. They use an auxiliary variable FLOW^® denoting an upper 
bound on the number of vehicles evacuated from area k. Constraints impose the 
capacity constraints. 


The model NEPP-SAT is the satisfaction problem version of NEPP-MF where the 
objective Q has been removed and the constraint 

FLOWfe = dk yk e £ (6) 

has been added to ensure that every vehicle is evacuated. 

To minimize clearance time, i.e., to find the minimal scheduling horizon such that 
all vehicles are evacuated, it suffices to add the objective function to NEPP-SAT 

min OBJECTIVE = max(END/c + tp^) (7) 

kGS 

and to relax the start and end time domains to [0, "H"®] where is an upper bound on 
the horizon required to evacuate all vehicles to a shelter. The resulting constraint-based 
scheduling model is denoted by NEPP-CT. 

A constraint-programming formulation NPEPP-MF of the non-preemptive phased 
evacuation planning problem can be obtained from NEPP-MF by replacing Q with ([T]|, 
which prevents the flows from two distinct origins to travel on the same edge at the same 
time. NPEPP-SAT, which is the satisfaction problem version of NPEPP-MF, is obtained 
by removing the objective Q and adding the constraint (|^. NPEPP-CT, which mini¬ 
mizes the evacuation clearance time, is obtained from NPEPP-SAT by adding the ob¬ 
jective 0. Note that since the flow-rate bounds ensure that edges capacities are always 
respected in NPEPP, the flow-rate variable can be directly set to its upper bound to max¬ 
imize evacuation efficiency. Hence, the following constraints are added to the NPEPP 
model: 

ApLowfe = inin((ife,Mp^) Vfc S £. (8) 

3.4 Problem Decomposition 

This section shows how the NEPP and NPEPP can be decomposed by understanding 
which paths compete for edges and/or how they relate to each other. In the following 
we introduce the path dependency relationship. 

Definition 1. Two paths px and Py are directly dependent, which is denoted by px^spy, 
if and only if they share at least a common edge, i.e., Ap„ , n Ap.^ 0. 

Definition 2. Two paths Px and Pz are indirectly dependent if and only if ^{pxApz) 
and there exists a sequence of directly dependent paths py .^,... such that px££py^, 

Py^Apy^ andPy^Apy.^,. . . ,Py^_^Apy^. 

Definition 3. Two paths px and Py are dependent, which is denoted by Px~^Py, if they 
are either directly dependent or indirectly dependent. Conversely, paths px and Py are 
independent, which is denoted by px-Lpy, if they are neither directly dependent nor 
indirectly dependent. 

Obviously, the path dependency T forms an equivalence relationship, i.e., T is reflex¬ 
ive, symmetric, and transitive. 



(a) Directly dependent paths 



(b) Indirectly dependent paths 


Fig. 2: Illustrating independent nodes and paths. 


The key idea behind the decomposition is to partition the evacuation areas E into 
T = {I? 0 ) ■ • ■ in such a way that any two paths, respectively from the set of 
evacuation areas Vi and Vj {0 < i < j < n), are independent. As a result, it is possible 
to solve the overall model by solving each subproblem Vi (0 < i < n) independently 
and concurrently. Figure illustrates two sets of evacuation nodes Vq = {x,y} and 
Vi = {z} where paths px and py are directly dependent and there are no indirectly 
dependent paths. Figure illustrates a single set of evacuation nodes V = {x,y,z} 
where the set of paths {px, Py} and {py, p^} are directly dependent, while the set of 
paths {px,Pz} are indirectly dependent. We now formalize these concepts. 

Definition 4. Let f?®* denote the paths of the set of nodes Vi. Two sets Vq and Vi 
of evacuation areas are independent if and only if any two paths from and 
respectively are independent. They are dependent otherwise. 

Theorem 1. Let T = {Vq, ..., Vn} be a partition of £ such that Di and Dj (0 < i < 
j < n) are independent. Then the NEPP and NPEPP can be solved by concatenating 
the solutions of their subproblems Vi (0 < i < n). 

The partition T — {Vq, ... ,!?„} can be generated by an algorithm computing the 
strongly connected components of a graph. Let Q be the directed graph consisting of 
the edges and vertices of all paths pk G and let be its undirected counterpart, 
i.e., the graph obtained after ignoring the direction of all edges in ■ The strongly 
connected components of define the partition T. 

3.5 Dominance Relationships 

This section shows how to exploit dominance relationships to simplify the constraint- 
based scheduling models. The key idea is to recognize that the capacity constraints of 
some edges are always guaranteed to be satisfied after introducing constraints on other 
particular edges. 

Definition 5. Let A a set of edges and e,e' ^ A, e e'. Edge e dominates d, denoted 
by e > e', if and only if 










- For simultaneous evacuation, the capacity of e is less than or equal to the capacity 

of e': < u^,; 

- The set of paths using e' is a subset of the set of paths using e: £{e') C £{e); 

- For non-convergent evacuation paths, the travel times for evacuation paths in £ (e') 
between e and e! are the same. 

Note that two edges may be dominating each other. For this reason and without loss 
of generality, this paper breaks ties arbitrarily (e.g., by selecting the edge closer to a 
safe node as the dominating edge). Note also that the capacity condition is ignored for 
phased evacuation. 

Theorem 2. Let the set of dominating edges in A We can safely substitute A^ to A 

in (0 in NEPP-MF such that the cumulative constraints are only stated for dominating 
edges. Similar results hold for NEPP-CT/SAT, and for the disjunctive constraints in 
NPEPP-MF/CT/SAT. 

3.6 Additional Constraints to a Real-World Evacuation Scheduling Prohlem 

The flexibility of the constraint-based evacuation scheduling approach allows to easily 
include many constraints appearing in real-world evacuation scheduling. For example, 
each flow rate variable domain may be restricted to a subset of values only, in order to 
account for the number of door-knocking resources available to schedule the evacua¬ 
tion d. Other real-world evacuation constraints may restrict departure times in order 
to wait for a more accurate prediction of an upcoming disaster or to restrict evacuation 
end times to ensure that the last vehicle has left a certain amount of time before the 
disaster strikes the evacuation area. 

3.7 Complexity of Phased Evacuations with Convergent Paths 

When using convergent paths, phased evacuations may be solved in polynomial time. 

Theorem 3. Model NPEPP-MF can be solved in polynomial time for convergent paths 
if all evacuation paths share the same latest completion time at the safe node. 

Proof (Sketch). Using the decomposition method and the dominance criterion, each 
subproblem with at least two evacuation paths includes exactly one dominating edge 
e which is part of every evacuation path. An optimal schedule can then be obtained in 
two steps. The first step builds a preemptive schedule by a sweep over the time starting 
from the minimal earliest start time of the tasks on e and ending before the shared 
completion time. For each point in time, it schedules one eligible task (if existing) 
with the largest flow rate (ties are broken arbitrarily) where a task is eligible if the 
point in time is not after its earliest start time on e and it has not been fully scheduled 
before that time. Note if a task is almost fully scheduled except the last evacuation 
batch then this step considers the actual flow rate of the last batch instead, which may 
be smaller than the task flow rate. This preemptive schedule is optimal as for each point 
in time, the unscheduled eligible tasks do not have an (actual) greater flow rate than 
the scheduled ones. The second step converts this optimal preemptive schedule to a 


non-preemptive one by postponing tasks interrupting others until after the interrupted 
tasks are completed. This transformation does not change the flows and hence the non- 
preemptive schedule has the same objective value as the preemptive schedule. □ 


3.8 The Search Procedure 

The search procedure considers each unassigned task in turn and assigns its under¬ 
lying variables. A task TASK^ is unassigned if the domain of any of its variables 
{STARTfc, DURfc, ENDfc, FLOWfc, ApLowfc} has more than one value. The search proce¬ 
dure selects an unassigned task TASK^ and then branch on all its underlying variables 
until they all are assigned. Depending on the considered problem, the models use differ¬ 
ent heuristics to (1) And the next unassigned task and to (2) select the next task variable 
to branch on and the value to assign. 

For NPEPP, the flow rate is directly set to the maximal value. The search strategy is 
determined by the problem objective as follows. If the objective maximizes the number 
of evacuees for a given scheduling horizon the search is divided into two steps. The 
first step selects a task with an unfixed duration and the highest remaining actual flow 
rate. If the lower bound on duration of the task is at least two time units less than its 
maximal duration then a minimal duration of the maximal duration minus 1 is imposed 
and a maximal duration of the maximal duration minus 2 on backtracking. Otherwise 
the search assigns duration in decreasing order starting with the largest value in its 
domain. The second step selects tasks according to their earliest start time and assigns 
a start time in increasing order|^ If the objective is to minimize the horizon % such 
that all vehicles are evacuated, then the search selects the next unassigned task with 
earliest start time by increasing order among all dominating edges, selecting the one 
with maximal flow rate to break ties, and to label the start time in increasing order. 

ForNEPP, the different search heuristics are as follow. Eor the choice (1), the strat¬ 
egy (lA) randomly selects an unassigned task, and performs geometric restarts when 
the number of backtracks equal to twice the number of variables in the model, using a 
growth factor of 1.5. The strategy (IB) consists in selecting the next unassigned task in 
decreasing order of evacuation demand for the dominating edge with the greatest num¬ 
ber of tasks. Eor the choice (2), the strategy (2A) first labels the flow rate in increasing 
order, then the task start time also in increasing order and, finally, the flow in decreasing 
order. The strategy (2B) first labels the flow rate in decreasing order, then the flow in 
decreasing order again and, finally, the start time in increasing order. 


4 Experimental Results 

This section reports experiments on a set of instances used in ||6l. These instances are 
derived from a real-world case study: the evacuation of the Hawkesbury-Nepean (HN) 
floodplain. The HN evacuation graph contains 80 evacuated nodes, 5 safe nodes, 184 
transit nodes, 580 edges and 38343 vehicles to evacuate. The experimental results also 

Note that the search procedure does not assume convergent paths or restrictions on the latest 
arrival times at the safe node. 



Table 1; The strongly connected components associated with each HN80-Ix instance. 
Instance #vehicles #scc see details 


HN80 38343 5 

HN80-I1.1 42183 4 

HN80-I1.2 46009 3 

HN80-I1.4 53677 5 

HN80-I1.7 65187 4 

HN80-I2.0 76686 4 

HN80-I2.5 95879 5 

HN80-I3.0 115029 5 


{22,9048}, {17,10169}, {14,6490}, {22,9534}, {5,3102} 
{1,751}, {2,1281}, {40,23656}, {37,16495} 

{2,1398}, {35,18057}, {43,26554} 

{2,1631}, {28,16737}, {27,19225}, {4,3824}, {19,12260} 
{22,16992}, {2,1980}, {42,33240}, {14,12975} 

{15,13974}, {38,40612}, {2,2330}, {25,19770} 

{32,36260}, {6,11214}, {16,17983}, {6,9324}, {20,21098} 
{5,11574}, {12,14184}, {19,23403}, {7,13068}, {29,39651} 


consider a class of instances HN80-Ix using the HN evacuation graph but a number of 
evacuees scaled by a factor x G {1-1,1.2,1.4,1.7, 2.0,2.5, 3.0} to model population 
growth. For simplicity, the experiments did not consider a flood scenario and assume 
that network edges are always available within the scheduling horizon H. It is easy to 
generalize the results to various flood scenarios. 

For each evacuation instance, a set of convergent evacuation paths was obtained 
from the TDFS approach 0 . The TDFS model is a MIP which is highly scalable as it 
aggregates edge capacities and abstracts time. The paths were obtained for the maxi¬ 
mization of the number of evacuees within a scheduling horizon of 10 hours, allowing 
preemptive evacuation scheduling. Thus, for each instance, the set of evacuation paths 
can be thought of as forming a forest where each evacuation tree is rooted at a safe 
node and each leaf is an evacuated node. In this particular case, each evacuation tree is a 
strongly connected component. It is important to emphasize that paths are not necessar¬ 
ily the same for the different HN instances, nor necessarily optimal for non-preemptive 
scheduling, which explains some non-monotonic behavior in the results. Table [T] re¬ 
ports the evacuation paths, the number of vehicles to evacuate (#vehicles), the number 
of strongly connected components (#scc), the number of evacuated nodes, the number 
of vehicles per see (see details) for each HN80-Ix instance. Each strongly connected 
component is represented by a pair {a;, y} where x is the number of evacuated nodes 
and y the number of vehicles to evacuate. 

The experimental results compare the flow scheduling results obtained with the 
NEPP and NPEPP approaches and the flow scheduling problem (PSP) formulation pre¬ 
sented in 0 . The PSP is solved using a LP and it relaxes the non-preemptive con¬ 
straints. Indeed, the flow leaving an evacuated node may be interrupted and restarted 
subsequently at any time t G H, possibly multiple times. Moreover, the flow rates in 
the PSP algorithm are not necessarily constant over time, giving substantial scheduling 
flexibility to the PSP but making it very difficult to implement in practice. Once again, 
the PSP comes in two versions. The objective of the core PSP is to maximize the number 
of vehicles reaching safety, while the objective of PSP-CT is to minimize the evacuation 
clearance time. In order to compare the PSP algorithm and the constraint-programming 
approaches of this paper fairly, the PSP is solved with a time discretization of 1 minute. 
The experiments for the NEPP and NPEPP models use different search heuristics and 
each experimental run was given a maximal runtime of 1800 seconds per strongly con- 





Table 2; Percentage of Vehicles Evacuated with ESP, NEPP-ME/S AT, NPEPP-ME/S AT. 




FSP 


NEPP-MF/SAT 

NPEPP-MF/SAT 

Instance 

CPU (s) 

Perc. 

Evac. 

CPU (s) 

Perc. 

Evac. 

Search 

CPU (s) 

Perc. 

Evac. 

HN80 

0.9 

100.0% 

0.2 

100.0% 

(SAT) 

{IB,2B} 

3.4 

96.9% 

HN80-I1.1 

1.1 

100.0% 

1538.9 

99.2% 

{lA, 2B} 

1295.9 

58.4% 

HN80-I1.2 

1.0 

100.0% 

0.4 

100.0% 

(SAT) 

{IB,2B} 

1444.4s 

57.7% 

HN80-I1.4 

1.3 

100.0% 

1347.5 

99.3% 

{lA, 2B} 

307.0 

73.0% 

HN80-I1.7 

1.8 

100.0% 

1374.9 

97.8% 

{lA, 2A} 

0.3 

59.0% 

HN80-I2.0 

2.0 

97.9% 

1770.1 

93.1% 

{lA, 2B} 

5.9 

52.8% 

HN80-I2.5 

1.8 

82.2% 

1664.1 

79.2% 

{lA, 2B} 

0.1 

51.5% 

HN80-I3.0 

1.4 

69.2% 

887.2 

67.5% 

{lA, 2B} 

0.1 

43.1% 


nected component. The results were obtained on 64-bits machines with 3.1GHz AMD 
6-Core Opteron 4334 and 64Gb of RAM and the scheduling algorithms were imple¬ 
mented using the programming language JAVA 8 and the constraint solver Choco 3.3.0, 
except for NPEPP-ME where the search was implemented in ObjectiveCP. 

Maximizing the Flow of Evacuees. Table [^compares, for each HN80-Ix instance and 
a 10-hour scheduling horizon, the percentage of vehicles evacuated (Perc. Evac.) and 
the solving time in seconds (CPU (s)) with PSP, NEPP-MP/SAT and NPEPP-MP/SAT. 
All solutions found with PSP are optimal and are thus an upper bound on the number 
of vehicles that can be evacuated with NEPP and NPEPP. Prior to solving NEPP-MP 
(resp. NPEPP-ME), the algorithm attempts to solve NEPP-SAT (resp. NPEPP-SAT) 
with a 60s time limit and, when this is successful, the annotation (SAT) is added next 
to the percentage of vehicles evacuated. As we make use of decomposition and parallel 
computing, the reported CPU for NEPP/NPEPP is the latest of the time at which the 
best solution is found among all strongly connected components. The table reports the 
best results across the heuristics, i.e., the run where the most vehicles are evacuated ; 
for the random strategy, the best result is reported across 10 runs (note that the stan¬ 
dard deviation for the objective value ranges between 0.4% and 1.1% only across all 
instances). The search strategy for the best run is shown in column (Search) as a com¬ 
bination {TaskVar, VarOrder} where TaskVar is the heuristic for choosing the next task 
variable and VarOrder is the heuristic for labeling the task variables. 

The results highlight the fact that the constraint-based simultaneous scheduling 
model finds very high-quality results. On the first five instances, with population growth 
up to 70%, the solutions of NEPP-MP are within 2.2% of the preemptive bound. This 
is also the case for the largest instance. In the worst case, the constraint-based schedul¬ 
ing model is about 4.9% away from the preemptive lower bound. It is thus reasonable 
to conclude that the constraint-based algorithms may be of significant value to emer¬ 
gency services as they produce realistic plans for large-scaled controlled evacuations. 
Por NPEPP-ME, the solver found optimal solutions and proved optimality for all in- 









Fig. 3; Quality of solutions over Time for NEPP-MF. 


stances, except HN80-I1.1 and HN80-I1.2 for which the best found solution was within 
0.1% of the optimal onej^The results indicate that a phased evacuation is much less ef¬ 
fective in practice and decreases the number of evacuees reaching safety by up to 40% 
in many instances. Unless phased evacuations allow an evacuation of all endangered 
people, they are unlikely to be applied in practice, even if they guarantee the absence of 
traffic merging. 

Figure shows how the quality of solutions improves over time for all HN-Ix in¬ 
stances which are not completely evacuated, for a particular run. For all instances, a 
high-quality solution is found within 10 seconds, which makes the algorithm applica¬ 
ble to a wide variety of situations. When practical, giving the algorithm more time may 
still produce significant benefits: For instance, on FIN-11.7 the percentage of vehicles 
increases from 93.0% to 97.6% when the algorithm is given 1800s. Such improvements 
are significant in practice since they may be the difference between life and death. 

Profile of the Evacuation Schedules. To demontrate the benefits of NEPP, it is useful 
to look at the evacuation profiles produced by the various algorithms. Recall that Pig.[2 
displays a repartition of departure times for seven evacuated nodes in the original HN80 
instance in the optimal solution produced by the PSP solver. The key observation is 
that, for several residential areas, the departure times are widely distributed within the 
scheduling horizon, indicating that the ESP solution makes heavy use of preemption. 
In the PSP solution, the number of vehicles departing at each time step is often equal 
to the path capacity. But there are also some suprising combinations {evacuated node, 
time step}, such as {3, 50}, {3, 84} and {3, 85} where the flow rate is respectively 
22, 3, and 12 for evacuation area 3. In summary, the PSP solution is unlikely to be the 
basis of a controlled evacuation: It is just too difficult to enforce such a complicated 
schedule. Figure shows a repartition of departure times for the same nodes in the 
original HN80 instance using the NEPP. The evacuation profile for the departure times 
is extremely simple and can easily be the basis of a controlled evacuation. Its simplicity 
contrasts with the complexity of the PSP solution and demonstrates the value of the 
constraint-programming approach promoted in this paper. 


^ In our experiments, the problem NPEPP satisfies the condition for Theorem]^ Thus, these in¬ 
stances can be solved almost instantly using the algorithm outlined in the proof of Theoremj^ 
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Fig. 4: Departure times and flows of 7 evacuated nodes for the NEPP solver. 


Table 3; Evacuation clearance time (CT) with FSP-CT, NEPP-CT, and NPEPP-CT. 


FSP-CT 


NEPP-CT 


NPEPP-CT 

Instance CPU (s) 

CT (min) 

CPU (s) 

CT (min) 

Search 

CPU (s) 

CT (min) 

HN80 4.3 

398 

1370.9 

409 

{IB,2A} 

0.2 

680 

HN80-I1.1 7.5 

582 

280.3 

616 

{lA, 2A} 

0.3 

1716 

HN80-I1.2 5.0 

577 

6.2 

590 

{lA, 2B} 

0.3 

1866 

HN80-I1.4 6.1 

587 

1386.0 

614 

{lA, 2A} 

0.4 

1226 

HN80-I1.7 7.3 

583 

1298.6 

610 

{lA, 2A} 

0.3 

2307 

HN80-I2.0 4.0 

625 

1713.7 

657 

{lA, 2B} 

0.4 

2909 

HN80-I2.5 8.4 

1092 

110.6 

1133 

{lA, 2B} 

0.4 

1884 

HN80-I3.0 9.4 

1232 

212.6 

1235 

{lA, 2B} 

0.3 

2467 


Minimizing the Clearance Time. Table[^compares, for each HN80-Ix instance, the min¬ 
imal clearance time in minutes (CT (min)) found with ESP-CT, NEPP-CT and NPEPP- 
CT. All solutions found with ESP-CT and NPEPP-CT are optimal for the given set of 
paths. Once again, solutions found by NEPP-CT are of high-quality and reasonably 
close to the preemptive lower bound produced by ESP-CT. In the worst case, the results 
are within 5.1% of the preemptive lower bound. The clearance times of the phased evac¬ 
uations, which are optimal, are significantly larger than for the NEPP. Note again that 
paths are different between instances and are not necessarily optimal with respect to 
different scheduling horizons, which explain inconsistencies such as the horizon found 
for HN80-I1.4 being shorter than the horizon found for HN80-I1.2 with NPEPP-CT. 


The Impact of the Flow Rates. The constraint-based scheduling models have the flow 
rates as decision variables, which increases the flexibility of the solutions. Tablej^stud- 
ies the benefits of this flexibility and compares the general results with the case where 
the flow rates must be selected from a specific set, here {2, 6,10,15, 20}. This is sim¬ 
ilar to the approach proposed in fTT), which uses a fixed set of response curves and 
their associated mobilization resources. Note that the column-generation algorithm in 
112 does not produce convergent plans and discretizes time. The results seem to indi¬ 
cate that flexible flow rates sometimes bring benefits, especially for the larger instances 
where the benefits can reach 3.0% ; nonetheless the possible loss when using fixed rates 
is not substantial and may capture some practical situations. 



















Table 4; Vehicles Evacuated with NEPP-MF with Flow Rates in {2,6,10,15, 20}. 


Instance 


CPU (s) 

Perc. Evac. 

Search 

HN80 


0.4 

100.0% (SAT) 

{IB, 2B} 

HN80-I1.1 


1538.9 

99.2% 

{lA, 2B} 

HN80-I1.2 


0.9 

100.0% (SAT) 

{IB, 2B} 

HN80-I1.4 


986.0 

99.5% 

{lA, 2B} 

HN80-I1.7 


1289.5 

97.1% 

{lA, 2A} 

HN80-I2.0 


1614.3 

91.0% 

{lA, 2B} 

HN80-I2.5 


1784.9 

77.0% 

{lA, 2B} 

HN80-I3.0 


1558.7 

65.6% 

{lA, 2B} 

Table 5; Comparison of ESP and NEPP problem sizes. 



FSP-10 

FSP-15 

NEPP-MF 

Instance 

#cols 

#rows 

#cols #rows 

#vars #ctrs 

HN80 

44651 

145880 

68651 218780 

1958 2288 


Comparison of Model Sizes. One of the benefits of the constraint-based scheduling 
models is that they do not discretize time and hence are highly scalable in memory re¬ 
quirements. This is important for large-scale evacuations which may be scheduled over 
multiple days. Table [^compares the ESP problem size for a scheduling horizon of 10 
hours (FSP-10) and 15 hours (FSP-15) with the NEPP-MF problem size for the HN80 
instance, when using 1 minute time steps. It reports the number of columns (#cols) and 
the number of rows (#rows) of the ESP model, as well as the number of variables (#vars) 
and the number of constraints (#ctrs) of the NEPP model. As can be seen, the number 
of variables and constraints grow quickly for the ESP model and are about 2 orders of 
magnitude larger than those in the NEPP-MF model which is time-independent. 

5 Conclusion 

This paper proposes, for the first time, several constraint-based models for controlled 
evacuations that produce practical and actionable evacuation schedules. These models 
address several limitations of existing methods, by ensuring non-preemptive schedul¬ 
ing and satisfying operational evacuation constraints over mobilization resources. The 
algorithms are scalable, involve no time discretization, and are capable of accommo¬ 
dating side constraints for specific disaster scenarios or operational evacuation modes. 
Moreover, the models have no restriction on the input set of evacuation paths, which can 
be convergent or not. Preliminary experiments show that high-quality solutions, with an 
objective value close to optimal preemptive solutions objectives, can be obtained within 
a few seconds, and improve over time. Future work will focus on improving the propa¬ 
gation strength of the cumulative constraint for variable durations, flows, and flow rates, 
and on generalizing the algorithm for the joint evacuation planning and scheduling. 
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